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Abstract 

This  paper  is  a  small  review  of  the  use  of  computer  simulations  and  especially  the  use  of  standard  quantum-mechanical  ab  initio  electronic 
structure  calculations  to  rationally  design  and  investigate  different  choices  of  chemicals/systems  for  lithium  battery  electrolytes.  Covered 
systems  and  strategies  to  enhance  the  performance  of  electrolytes  will  range  from  assisting  the  interpretation  of  vibrational  spectroscopy 
experiments  over  development  of  potentials  for  molecular  dynamics  simulations,  to  the  design  of  new  lithium  salts  and  the  lithium  ion 
coordination  in  liquid,  polymer,  and  gel  polymer  electrolytes.  Examples  of  studied  properties  include  the  vibrational  spectra  of  anions  and 
ion  pairs  to  characterize  the  nature  and  extent  of  the  interactions  present,  the  lithium  ion  affinities  of  anions,  important  for  the  salt  solvation 
and  the  ability  to  provide  a  high  concentration  of  charge  carriers,  the  HOMO  energies  of  the  anions  to  estimate  the  stability  versus  oxidation, 
the  anion  volumes  that  correlate  to  the  anion  mobility,  the  lithium  ion  coordination  and  dynamics  to  reveal  the  limiting  steps  of  lithium  ion 
transport,  etc. 
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1.  Introduction 

The  lithium  ion  containing  electrolyte,  whether  liquid, 
polymer,  or  gel  polymer  based,  is  perhaps  the  part  of 
the  lithium  batteries  of  today  where  there  is  most  room 
for  improvement.  Properties  of  the  bulk  electrolyte  may 
limit  the  battery  performance  (e.g.  degradation,  flammability, 
and  voltage  window),  as  may  the  properties  of  the  elec¬ 
trolyte/electrode  interfaces.  Especially,  the  nature  of  the  solid 
electrolyte  interfaces  (SEIs),  both  at  anodes  and  cathodes,  are 
strongly  dependent  on  the  electrolyte  composition.  There¬ 
fore,  there  is  a  strong  interest  to  develop  better  battery  elec¬ 
trolytes  and  preferably  through  a  fundamental  understanding 
of  the  factors  finally  limiting  battery  performance. 

However,  direct  structure-property  relations  are  trouble¬ 
some  to  get  about  for  electrolytes.  This  is  due  both  to  the 
range  of  possible  interactions  within  the  materials  and  in 
the  electrode/electrolyte  interfaces  as  well  as  the  disordered 
structure,  which  limits  the  applicability  of  experimental  struc¬ 
ture  characterization  techniques  severely  compared  to  those 
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available  for  studying  electrode  materials.  Simultaneously 
there  is  a  driving  force  to  use  more  complex  and  specialized 
electrolytes,  which  further  complicates  analysis  of  the  inter¬ 
actions  present  and  the  linkage  between  each  component  and 
the  resulting  electrolyte/battery  performance. 

There  are  also  high  costs  associated  with  the  development 
of  new  electrolytes.  It  is  thus  highly  essential  to  develop  a 
cost-effective  and  systematic  method  for  selection  of  appro¬ 
priate  electrolyte  components.  Ideally,  such  a  method  should 
be  used  as  a  tool  for  rational  design  of  new  components 
prior  to  usage  or  synthesis  efforts,  thereby  serving  as  an 
effective  and  reliable  screening  device.  It  should  also  be 
advantageous  with  an  as  general  and  easily  accessible  method 
as  possible  to  make  comparisons  using  data  from  different 
sources,  e.g.  research  groups  and  equipments,  unambiguous. 
With  the  advent  of  many  academic  as  well  as  commercial 
software  packages,  efficient  and  diversely  tailored  compu¬ 
tational  algorithms,  and  the  everywhere  abundant  computa¬ 
tional  resources  it  is  perhaps  natural  to  seek  the  answer  within 
publicly  available  standard  ab  initio  computational  chemistry 
methods. 

Computer  simulations  can  in  many  cases  be  successfully 
applied  to  electrolyte  materials.  The  purpose  and  outcome 
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of  the  simulations  differ  with  the  method  applied:  molecular 
dynamics  (MD)  and  Monte-Carlo  (MC)  simulations  model 
properties  accessible  using  large  ensembles.  MD  simulations 
have  been  applied  by  many  different  research  groups  to 
a  wide  range  of  non-aqueous  liquid  electrolytes  [1-5] 
and  to  the  more  complex  and  intriguing  solid  polymer 
electrolytes  (SPEs)  [6-25],  with  a  strong  emphasis  on 
poly(ethylene  oxide)  (PEO)-based  lithium  ion  containing 
ones  [6-11],  but  also  SPEs  based  on  other  polymers  than 
simple  linear  PEO  [10,12,13],  as  well  as  SPEs  based  on 
other  cations  than  lithium  ions  [8,14-25],  especially  sodium 
ions  [14-17,20-22,24,25].  Recently  also  crystalline  SPEs 
[26]  have  been  subject  to  MD  simulations  [27]  and  also 
polymer  electrolytes  containing  nano-particles  [28-30]. 
MC  simulations  have  been  less  frequently  used,  but  have 
revealed  the  limiting  step  in  ion  conduction  in  polymer 
matrices  by  the  development  and  testing  of  the  dynamic 
bond  percolation  theory  [31],  have  been  used  to  model  the 
Li+  adsorption  on  a  metal  electrode  from  a  liquid  electrolyte 
[32],  and  explained  the  difference  between  the  conduction 
mechanisms  in  polymer  electrolytes  and  polyelectrolytes 
[33,34]. 

On  the  other  hand,  ab  initio  calculations  can  only  be 
applied  to  much  smaller  systems  and  are  thus  directed  to 
probe  properties  dependent  on  local  phenomena.  The  major 
benefit  of  ab  initio  calculations  compared  to  MD  or  MC, 
apart  from  the  inclusion  of  electrons  and  thereby  electronic 
properties  directly,  is  that  the  results  are  independent  of  any 
previously  developed  force-field  and  that  local  phenomena 
are  treated  undisturbed.  In  this  paper,  we  will  review  the 
development  of  the  usage  of  ab  initio  calculations  with 
respect  to  lithium  battery  electrolytes.  We  have  identified 
four  different  areas/systems  where  ab  initio  calculations 
have  been  applied  to  study  properties  of  battery  electrolytes 
either  alone,  in  conjunction  with  experimental  techniques, 
or  to  serve  further  simulations  with  necessary  input  data. 
The  areas:  (1)  anions  and  ion  pairs  calculating  the  lithium 
ion  affinities  of  the  anions  and  the  effect  of  formation  of 
ion  pairs  in  vibrational  spectra  used  for  characterization  of 
electrolytes.  For  the  anions  also  the  volume  and  oxidative 
stability,  important  for  the  diffusion  and  the  electrochemical 
stability  window,  respectively,  can  be  calculated,  (2)  polymer 
conformation  and  flexibility — in  polymer  electrolytes  the 
local  structures  and  changes  within  the  polymer  chains  and 
the  associated  energies  and  energy  barriers  are  very  impor¬ 
tant  for  the  ion  conduction  mechanism,  (3)  the  lithium  ion 
coordination — the  nature  of  the  first  solvation  shell  of  solvent 
molecules,  polymer  chains,  or  combinations  of  these  around 
a  cation  core,  crucial  to  gain  a  molecular  level  insight  to  the 
cation  conduction  mechanism.  This  also  includes  the  lithium 
ion  environment  in  the  new  crystalline  polymer  electrolytes, 
LiYF6P(EO)6  (Y  =  P,  As,  and  Sb)  made  by  the  St.  Andrews 
group  [26],  and  (4)  additives  to  electrolytes  to  enhance 
the  safety  and/or  the  life-time  of  the  batteries,  often  by 
promoting/prohibiting  reactions  at  the  electrolyte/electrode 
interfaces. 


In  general,  the  use  of  ab  initio  methods  also  allows  the  col¬ 
lection  of  different  basic  physical  properties:  local  and  global 
minima  on  potential  energy  surfaces,  bond  strengths,  interac¬ 
tion  strengths,  charge  distributions,  etc.  Below  we  summarize 
the  data  and  interpretations  collected  over  the  years,  together 
with  an  outlook  on  future  possibilities  in  the  area  of  applying 
ab  initio  methods  to  battery  electrolytes. 

2.  Case  studies 

2.7.  Anions  and  lithium  ion-anion  pairs 

Ions  and  ion  pairs  were  early  studied  using  ab  initio  meth¬ 
ods,  mainly  for  simple  ions  like  the  halogens,  but  in  the  1990s 
work  on  larger  polyatomic  anions  and  their  lithium  ion  pairs 
started  to  appear,  among  the  first  being  a  study  on  LiBF4  [35]. 
In  1993,  the  first  paper  using  ab  initio  calculations  with  the 
outspoken  purpose  of  assisting  in  studies  of  lithium  battery 
electrolytes  was  published  on  the  triflate  (Tf),  (CF3SO3-), 
anion  [36],  quickly  followed  by  two  other  papers  on  the  tri¬ 
flic  acid  [37],  and  the  lithium  ion  pairs  [38].  The  following 
year  another  study  on  the  triflate  anion  also  including  lithium 
ion  pairs  and  higher  aggregates  appeared  [39]. 

These  studies,  on  the  at  the  time  very  popular  electrolyte 
anion  Tf,  all  focused  on  assisting  in  interpretation  of  vibra¬ 
tional  spectra  (IR  and  Raman).  Thus,  Gejji  et  al.  calculated 
the  vibrational  spectra  for  Tf,  for  its  corresponding  acid  HTf, 
as  well  as  for  mono-  and  bi-dentate  1 : 1  LiTf  ion  pairs.  In  the 
work  by  Huang  et  al.,  the  coordination  environment  of  the 
lithium  cation  was  altered  to  mimic  a  solvent  by  including 
OH-groups,  for  a  more  realistic  interaction  with  the  anion. 
One  controversy  among  experimentalist  was  the  order  of 
the  S-0  and  C-F  asymmetric  stretching  vibrations.  Together 
with  the  experimental  and  computational  work  on  180  sub¬ 
stituted  anions  by  Johnston  and  Shriver  [40],  the  assignment 
of  the  vibrational  spectra  of  the  “free”  Tf  anion  and  the  ion 
pairs  was  clarified,  making  future  spectroscopic  measure¬ 
ments  on  ionic  interactions  in  electrolytes  straight-forward 
for  this  anion/lithium  salt. 

Similar  controversies,  this  time  concerning  the  nature  of 
the  lithium  ion  pairs:  mono-,  bi-,  or  tridentate,  existed  for  the 
C104"  anion.  In  1996,  Klassen  et  al.  [41]  presented  a  paper 
clearly  directed  to  the  electrolyte  community  showing  that 
the  preferred  coordination  is  bi-dentate.  However,  Ramondo 
et  al.  had  made  the  same  conclusion  earlier,  albeit  using  lower 
level  calculations  [42] . 

Due  to  the  increasing  amount  of  computational  resources 
available,  also  larger  anions  become  possible  to  study,  fol¬ 
lowing  the  experimental  development  of  anions  with  more 
delocalized  charges.  Also  a  more  elaborate  study  on  Tf  was 
made  possible  [43].  In  1995,  a  comparative  study  on  the 
Tf,  TFSI,  [N(CF3S02)2_]  (a.k.a.  the  “imide”  anion),  and 
TFSM,  [CH(CF3S02)2_]  (a.k.a.  the  “methide”  anion),  com¬ 
mon  being  the  electron  withdrawing  CF3S02  group,  was 
published  [44].  This  study  was  directed  towards  the  charge 
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distribution,  analysis  of  the  binding  within  the  anions  in 
detail,  and  to  reveal  stabilities  by  computing  the  chemical 
hardness  from  the  HOMO  and  LUMO  energies. 

Later,  a  potential  energy  surface  (PES)  investigation 
revealed  two  possible  stable  conformers  of  the  TFSI  and  its 
large  internal  flexibility,  providing  an  additional  explanation 
for  the  observed  plasticizing  effect  apart  from  its  delocal¬ 
ized  charge  [45].  Also,  the  vibrational  spectra  for  the  TFSI 
anion  and  the  HTFSI  acid  was  computed,  revealing  earlier 
errors  in  the  assignment  of  bands  used  for  observation  of  ion 
pairing  in  the  electrolytes  [46] .  Finally  also  the  ion  pairs  of 
FiTFSI  and  FiTFSM  have  been  computed  [47,48],  reveal¬ 
ing  the  common  nature  of  the  binding  to  two  S02-groups  for 
TFSI  and  TFSM,  but  showing  also  that  TFSI  can  complex  the 
lithium  cation  by  the  central  nitrogen  atom  [48].  A  PES  study, 
very  similar  to  the  TFSI  one,  was  also  made  for  the  PFSI  (or 
“beti”)  anion  [49].  In  the  late  1990s  and  early  2000s,  it  was 
possible  to  compute  a  wealth  of  properties  for  almost  any 
anion/ion  pair  of  interest  at  a  reasonable  computational  level. 
Example  of  such  studies  using  large  anions  are  the  TFSM 
and  CTFSM  [C(CF3S02)3~]  anions  and  their  ion  pairs  [50], 
aromatic  lithium  salts  [51],  and  the  series  of  salts  with  the 
formula  FiPF6-n(CF3)n  [52]. 

We  identify  today  two  development  lines  of  applying  ab 
initio  methods  to  anions  and  ion  pairs:  re-calculating  earlier 
systems  with  much  higher  accuracy  or  extending  the  systems 
to  include  also  solvents  by  explicit  molecules  or  by  using 
continuum  models.  Examples  of  the  former  are  the  high  accu¬ 
racy  calculations,  using  very  large  and  accurate  basis  sets, 
to  develop  better  Fi+-PF6_  potentials  for  MD  simulations 
[53],  studies  of  higher  aggregates  of  [Fix+(AsF6_)>;]  [54],  a 
new  PES  study  on  TFSI  (and  TFSM)  [55]  and  using  DFT 
methods  to  correlate  dissociation  and  lattice  energies  [56]. 
Examples  of  the  latter  are  the  study  on  FiBF4  in  acetoni¬ 
trile  using  explicit  solvent  molecules  [57],  the  work  by  Wang 
and  Balbuena  on  FiTf  in  polyphosphazene  and  water  [12], 
the  study  on  crystalline  solvates  using  FiPF6  and  FiCICL 
and  ethers  by  Grondin  et  al.  [58],  the  of  different  ion  pairs 
together  with  a  solvent  molecule  by  us  [59] ,  and  studies  on  ion 
pairs  using  self-consistent  reaction  field  continuum  methods 
[2,60].  Also,  calculations  to  support  interpretation  of  vibra¬ 
tional  spectra  as  before  continues,  but  now  applied  to,  e.g.  the 
new  salts  FiTADC  [Fi+(N5C4“)]  andFiPATC  [Fi+(N5C6_)] 
based  on  heterocyclic  anions  [61-63]. 

However,  we  find  the  most  exciting  development  of  using 
ab  initio  methods  the  possibility  to  predict  not  only  sta¬ 
ble  structures  and  their  vibrational  spectra  to  assist  in  the 
characterization  of  electrolytes,  but  rather  to  assist  in  the 
development  of  new  salts  by  assessing  the  (relative)  inter¬ 
action  strengths.  This  strategy  is  inspired  by  Kim  et  al., 
showing  the  excellent  correlation  between  salt  lattice  energies 
and  interaction  energies  for  1:1  cation-anion  systems  [56]. 
Hereby,  anions  with  low  lithium  ion  affinities  can  be  devel¬ 
oped  “in  silico”  (by  use  of  computations  prior  to  synthesis 
work)  in  a  reliable  manner,  saving  time  and  money.  In  Fig.  1, 
the  obtained  electronic  energy  differences  for  the  reaction: 
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Fig.  1.  Electronic  energy  differences  between  free  anions  and  cations  vs. 
ion  pairs  for  families  of  lithium  salts. 


FiAn  Fi+  +  An- ,  for  different  classes  of  anions  are  plotted. 
The  anions  are  divided  into  seven  families,  albeit  somewhat 
ambiguously:  “simple”  anions  are,  e.g.  mono-atomic  and 
OCN“,  HC02“,  N03“,  SCN“,  CH3COO-,  etc.,  “classic” 
weakly  coordinating  anions  are,  e.g.  BF4-,  PF^-,  CIO4-, 
Tf,  etc.,  “imides”  are  TFSI,  TFSM,  PFSI,  etc.,  “B-based” 
anions  those  with  a  boron  atom  center,  e.g.  bis-(oxalate)- 
borate  (BOB),  malonato-oxalato-borate  (MOB),  etc.,  “P- 
based”  anions  those  with  a  phosphorous  atom  center,  e.g. 
PF6_„(CF3),T,  PF3(C2F5)3-  (“FAP”),  etc.,  “heterocyclic”, 
e.g.  TADC,  PATC,  etc.,  and  “other”,  e.g.  the  N(CN)2-  anion. 
These  are  original  data  based  on  new  calculations,  details  of 
which  are  to  be  presented  elsewhere  [64],  though  using  start¬ 
ing  geometries  from  the  literature  to  a  large  extent.  On  the 
basis  of  these  data,  there  are  possibilities  to  focus  on  the  types 
of  anion  chemistries  most  promising  with  respect  to  reduc¬ 
ing  the  lithium  ion  affinity,  eventually  resulting  in  electrolytes 
with  higher  charge  carrier  concentrations. 

The  calculations  can  also  be  used  to  evaluate  the  stabil¬ 
ity  versus  oxidation  for  each  anion  using  the  HOMO  energy 
(Fig.  2)  as  previously  described  [65,66].  The  volume  of  the 
anion  is  another  important  attainable  parameter  as  the  anion 
diffusion  in  the  electrolyte  matrix  should  be  inversely  pro- 


Fig.  2.  HOMO  energies  for  families  of  anions. 
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Fig.  3.  Volumes  for  families  of  anions. 


portional  to  the  radius  of  the  anion.  A  slow  diffusion  due  to 
a  large  anion  should  decrease  the  total  ion  conductivity,  but 
increase  the  cation  transference  number — one  of  the  most 
crucial  goals  to  resolve  for  lithium  battery  electrolytes.  The 
extreme  is  the  anion  tethered  or  incorporated  into  the  matrix 
as  in  single  cation  conducting  poly  electrolytes.  In  Fig.  3,  the 
volumes  of  the  different  families  of  anions  are  plotted  using 
the  electronic  spatial  extent  as  volume  measure.  The  wide 
range  of  volumes  available  for  most  anion  families  should 
allow  to  specifically  target  the  balance  total  ion  conductivity 
and  cation  transference  numbers  a  priori. 

By  studying  the  data  further,  we  find  the  naive  notion  of 
a  correlation  between  a  large  anion  and  a  highly  dissociative 
salt  (Fig.  4),  to  even  qualitatively  only  be  approximate  for 
very  small  anions/large  electronic  energy  differences. 

Together  these  data  provide  information  to  justify  choices 
of  anions/lithium  salts  depending  on  preferences  for  the  bat¬ 
tery  electrolyte,  e.g.  voltage  and  temperature  working  range, 
cation/anion  transport  number  etc,  hence  showing  the  power 
of  ab  initio  calculations  on  anions  and  ion  pairs  to  assist  in 
rational  design  of  electrolytes. 

One  practical  way  is  to  apply  quantitative  structure- 
property  relationships  (QSPR)  to  the  data.  The  choices  of 


Fig.  4.  Volumes  of  anions  vs.  electronic  energy  differences. 


independent  targets  are,  e.g.  small  electronic  energy  differ¬ 
ences  (as  defined  above  for  Fig.  1),  low  HOMO  energies 
(large  negative  values);  hence  aiming  at  highly  dissociative 
lithium  salts  that  are  electrochemically  stable.  Anion  descrip¬ 
tors  can  be  based  both  on  primary  chemical  structure  (#atoms, 
#halogens,  #0,  #N,  etc.)  and  on  computed  properties  (anion 
size,  HOMO,  LUMO,  maximum  atomic  charge,  etc.).  Exam¬ 
ples  of  ion  pair  descriptors  are:  electronic  energy  differences, 
Li+  coordination  number,  etc.  A  tentative  QSPR  analysis 
shows  highly  dissociative  salts  to  be  products  of  a  large  anion, 
a  low  LUMO  energy  for  the  anion,  and  preferably  a  Li+-N 
coordination  site.  Low  HOMO  energies,  characterizing  elec¬ 
trochemically  stable  anions,  are  by  the  QSPR  analysis  found 
to  be  products  of  high  LUMO  energies  and  many  halogen 
atoms  in  the  structure. 

However,  the  limited  range  of  anions  heavily  affects  the 
results  and  thus  these  conclusions  should  be  used  with  care. 
Though,  it  is  re-assuring  that  the  QSPR  results,  e.g.  that 
heavily  fluorinated  anions  are  electrochemically  more  stable 
versus  oxidation  and  large  anions  in  general  provide  disso¬ 
ciative  salts,  agree  with  knowledge  from  experimental  data. 

2.2.  Polymer  conformation  and  dynamics 

While  organic  solvents  for  liquid  electrolytes  easily  can  be 
studied  by  ab  initio  calculations,  polymers  for  solid  polymer 
electrolytes  are  special  and  have  properties  difficult  to  grasp 
experimentally  and  thus  merit  an  ab  initio  treatment.  Simulta¬ 
neously,  polymers  are  difficult  to  attack  by  ab  initio  methods, 
making  the  scientific  challenge  greater.  Amongst  the  poly¬ 
mers  used  poly  (ethylene  oxide),  PEO  (-CH2CH2O-L,  has 
been  the  archetype  ever  since  the  first  paper  on  alkali  metal 
polymeric  electrolytes  [67]  appeared  in  1973.  Therefore,  the 
dominant  model  systems  used  in  ab  initio  calculations  to 
model  SPE  polymer  conformation  have  been  oligomers  of 
PEO  and  often  the  so-called  glycol  dimethyl  ethers  (“gly- 
mes”),  CH30(CH2CH20)nCH3.  The  use  of  a  methyl  end- 
capped  oligomer  is  a  very  good  approximation  to  high  Mw 
PEO  {n  ~  105)  with  two  OH-endgroups.  Although  oligomer 
models  cannot  be  expected  to  reveal  the  ion  conduction 
behaviour  of  a  polymeric  system,  they  can  still  be  used  to 
investigate  the  local  dynamics  and  energies  related  to  con¬ 
formational  changes. 

Monoglyme  (CH30(CH2CH20)iCH3)  has  been  exten¬ 
sively  used  by  several  groups  for  studies  of  PEO  confor¬ 
mations  [68-80].  Its  small  size  allowed  for  treatment  with 
high  accuracy  models  (MP2  and  higher)  and  large  basis  sets 
already  in  the  1980s,  and  use  of  continuum  environments 
in  the  1990s  [68].  The  main  drawback  of  monoglyme  is  the 
short  length  making  no  part  of  the  oligomer  truly  “non-end”. 
Several  studies  have  specifically  been  aimed  at  calculating 
preferred  conformations  and  internal  torsion  barriers  between 
these  [68-73,75,79-81].  The  agreement  between  the  size  of 
the  energy  barriers  in  PEO  with  internal  rotations  in  TFSI 
[45],  revealed  the  specific  plasticizing  effect  of  the  latter  in 
LiTFSIPEOft  electrolytes. 
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Only  one  similar  investigation  for  the  popular  experi¬ 
mental  model  polymer  system  poly(propylene  oxide)  (PPO) 
exists  to  the  best  of  our  knowledge  [82].  This  may  be  due 
to  the  fact  that  the  “PPO-monoglyme”  has  a  higher  compu¬ 
tational  cost  and  an  increased  amount  of  conformers  due  to 
the  tactic/atactic  possibilities.  The  work  by  Sasanuma,  how¬ 
ever,  is  an  extensive  and  clear  account  with  most  aspects  of 
4  ‘PPO-monoglyme’  ’ . 

One  topic  of  interest  from  the  polymer  electrolyte  com¬ 
munity,  as  it  may  also  effect  the  salt  solvation  capacity,  has 
been  the  question  if  ab  initio  methods  can  explain  the  so- 
called  gauche  effect  of  O-C-C-O-based  polymers,  i.e.  why 
the  gauche  conformation  about  a  C-C  bond  is  as  or  even 
more  stable  than  an  anti  conformation  [68,74,77,80,82].  The 
effect  has  been  verified  both  by  using  infrared  measurements 
in  argon-matrices  [69]  and  electron  diffraction  [83],  and  one 
suggestion  made  is  the  presence  of  a  non-bonded  C-H-  •  0 
preferential  contact,  but  there  is  still  no  consensus  as  to  the 
origin  of  the  effect  nor  to  the  resulting  properties. 

The  next  oligomer  in  length  -  the  diglyme,  CH30(CH2- 
CH20)2CH3  -  has  also  found  extensive  use,  although  its 
larger  size  for  a  long  time  prohibited  use  of  more  sophisticated 
methods  than  Hartree-Fock.  In  1994,  Gejji  et  al.  showed  the 
effect  of  end-group  versus  non-end-group  to  be  crucial  when 
generating  torsion  potentials  [84].  Diglyme,  having  both 
O-C-C-O  and  C-O-C-C/C-C-O-C  dihedral  angels  that 
do  not  involve  atoms  inherently  being  from  the  methyl  end- 
groups,  was  found  to  be  a  much  better  PEO  model  than  mono- 
glyme.  In  1997,  a  thorough  investigation  with  large  basis  sets 
and  MP2  level  calculations  was  performed,  which  revealed 
a  large  number  of  possible  stable  conformations  [85].  The 
barriers  for  conformational  changes  were  re-calculated  and 
the  authors  suggested  the  gauche  effect  to  become  more  pro¬ 
nounced  for  longer  oligomer  chains.  Torsion  potentials,  for 
MD  simulations,  have  in  many  cases  been  obtained  by  fitting 
analytical  expressions  to  values  from  single-point  ab  initio 
calculations  [15-17].  This  is  one  way  to  get  high  accuracy 
MD  potentials  for  a  specific  system,  which  for  polymers-like 
PEO  might  be  crucial  as  the  intermolecular  conversions,  i.e. 
conformational  changes,  have  low  energy  barriers. 

Following  a  decreased  interest  in  PEO-based  electrolytes 
and  other  SPEs,  together  with  the  wealth  of  data  already 
accessible,  the  activity  within  ab  initio  calculations  to  reveal 
polymer  flexibility  has  declined,  though  new  MD  potentials 
still  are  published  [86]. 

2.3.  Lithium  ion  coordination 

Whether  the  environment  is  liquid,  polymeric,  gel-like,  or 
solid,  understanding  the  nature  of  the  lithium  ion  coordination 
is  crucial  to  understanding  the  conductivity  mechanism.  All 
the  above  situations  have  to  different  extents  been  attacked 
by  ab  initio  calculations. 

For  lithium  ion  batteries  and  liquid  electrolytes,  cyclic 
organic  solvents  like  tetrahydrofuran  (THF),  ethylene  car¬ 
bonate  (EC),  propylene  carbonate  (PC),  y-butyrolactone 


(GBL),  and  their  mixtures  have  been  proven  to  be  efficient, 
especially  with  respect  to  cyclability.  A  mixture  of  two  or 
more  solvents  allows  the  design  of  different  electrolyte  prop¬ 
erties  in  a  wider  range  and  thus  to  enhance  performance.  The 
application  of  ab  initio  calculations  to  liquid  electrolytes  for 
lithium  batteries  was  pioneered  by  Blint  and  his  calculations 
on  different  ether  and  carbonyl  oxygen  containing  species  and 
their  mixtures  and  the  coordination  of  lithium  ions  [87,88]. 
His  results  on  EC/PC  coordination  of  Li+,  always  preferen¬ 
tially  by  the  carbonyl  oxygen  atoms,  have  later  been  both 
verified  and  further  refined  by  several  authors  [2,4,89-91]. 
Blint  extrapolated,  using  data  on  three  EC  ligands  around 
Li+  that  a  four-coordinated  complex  would  be  dominant. 
Klassen  et  al.  computed  exothermal  AH  for  complexes  with 
up  to  four  EC  ligands  and  recent  higher  levels  of  calculations 
by  Wang  et  al.  showed  the  Gibbs  free  energy  to  support  the 
[Li(EC)4]+  complex,  with  ^S4  symmetry,  to  be  the  prevail¬ 
ing  component.  The  data  are  in  excellent  agreement  with  MD 
simulations  [1,2]  and  spectroscopic  investigations  [92,93], 
which  further  showed  no  preferential  solvation  in  EC/PC 
mixtures,  making  the  ab  initio  results  using  only  EC  ligands 
valid  also  for  EC/PC  mixtures.  Apart  from  the  studies  by 
Blint,  which  considered  also  di-methyl  ether,  di-ethyl  ether, 
acetone  and  water  [87],  and  acetaldehyde  [88],  some  recent 
studies  have  used  acetonitrile  [57]  and  y-butyrolactone  [5]  as 
solvents.  The  latter  solvent  has  gained  interest  due  to  its  wide 
liquid-phase  temperature  range  (—42  to  206  °C). 

For  modeling  the  lithium  ion  coordination  in  polymer  sol¬ 
vents  the  preference  for  PEO  model  compounds  is  clear. 
Several  ab  initio  studies  on  1:1  complexes  of  monoglyme 
and  a  lithium  ion  were  made  in  the  mid  1990s  [6,94-98].  The 
results  were  correlated  with  experimental  studies  on  lithium 
salts  both  in  PEO  and  in  monoglyme  [99,100].  Following 
the  computer  development  larger  model  compounds  for  PEO 
have  been  used:  diglyme  (in  both  1 : 1  and  1:2  complexes  with 
a  lithium  ion)  [101-105],  triglyme  [106],  and  longer  glymes 
[107,108].  First  with  the  use  of  longer  glyme  models  a  maxi¬ 
mum  coordination  number  of  six  for  the  lithium  cation  could 
be  obtained.  On  the  other  hand,  the  di-  and  tri-glyme  model¬ 
ing  efforts  were  successful  in  discerning  whether  the  lithium 
ion  transport  mechanism  in  a  polymeric  matrix  involves  the 
anion  or  not  [104],  whether  ion  conduction  is  due  to  SnI- 
or  SN2-type  mechanisms  [102-104],  and  located  transition 
states  of  lithium  ion  movement  along  polymer  chains  and 
the  associated  barriers  along  the  dominant  internal  reaction 
coordinate  paths  [102-104].  It  seems  as  if  higher  Li+  coor¬ 
dination  numbers  reduces  the  migration  barriers  regardless 
whether  or  not  an  anion  is  present  [104]. 

Recently,  other  polymer  matrices  have  been  studied, 
among  them  the  PEO-like  poly  alkyl  oxides,  PPO  and  PTMO, 
poly(trimethylene  oxide)  [109,110],  and  polyphosphazene 
(PP)  [12].  The  two  former  studies  revealed  steric  hindrances 
to  cause  the  stronger  complexation  by  PEO  compared  to  the 
other  polymers  for  a  six-coordinated  lithium  cation,  while  the 
latter  showed  the  importance  on  the  PP  backbone  nitrogen  in 
the  cation  complexation.  If  other  polymers  are  exploited  as 
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matrices,  there  are  today  no  obstacles  from  the  computational 
side  to  study  the  local  coordination  of  the  cation,  regardless 
of  the  size  of  the  monomer  unit. 

A  totally  different  polymeric  environment  for  the  lithium 
cation  is  found  in  the  crystalline  polymer  electrolyte 
LiYF6P(EO)6  (Y  =  P,  As,  and  Sb)  [26].  Recently,  we  inves¬ 
tigated  the  lithium  ion  coordination  and  path  of  migration 
within  LiPF6P(EO)6  by  using  single-point  energy  calcula¬ 
tions  on  an  oligomer  model  [111].  The  obtained  energy  bar¬ 
rier,  1  eV,  compares  favorably  with  the  Arrhenius  activation 
energy  obtained  from  ion  conductivity  data  (~1  eV). 

Gel  electrolytes  have  been  commercially  used  since  1999 
and  are  today  the  dominating  lithium  battery  concept  [112]. 
In  gels,  the  macroscopic  properties  of  a  polymer  is  com¬ 
bined  with  the  liquid-like  solvent  properties,  and  this  may 
affect  also  the  lithium  ion  coordination  environment.  Exam¬ 
ples  of  polymers  often  used  are  poly(acrylonitrile)  (PAN) 
and  poly  (methyl  methacrylate)  (PMMA).  Examples  of  sol¬ 
vents  used  are  ethylene  carbonate,  propylene  carbonate,  or 
mixtures  thereof.  Whether  the  PAN  or  PMMA  polymers 
take  an  active  part  in  coordinating  the  lithium  cation  can 
be  studied  by  computing  binding  energies  using  different 
combinations  or  polymer  and  solvent  molecules  [91].  For 
a  PMMA-based  gel,  the  pure  solvent  complex,  [Li(EC)4]+, 
seems  totally  dominant,  in  excellent  agreement  with  dis¬ 
cussions  based  on  spectroscopic  data  where  the  interaction 
with  PMMA  is  suggested  to  be  very  weak  [91-93].  For  a 
PAN-based  gel  a  significant  Li+-polymer  interaction,  which 
increases  with  temperature,  was  found  both  computationally 
and  spectroscopically  (Fig.  5)  [91]. 

A  comparison  not  made  before  is  that  of  cation  coordina¬ 
tion  across  liquid  (e.g.  EC/PC)  to  gel  (e.g.  PAN/EC/PC)  to 
polymer  (e.g.  PEO)  systems.  We  find  that  all  three  systems 
provide  about  equal  binding  energies  for  the  lithium  cation, 
despite  the  different  maximum  coordination  numbers  of  Li+ 
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Fig.  5.  Raman  spectra  of  the  LiC104/EC/PC/PAN  gel  polymer  electrolyte 
in  the  C-N  stretching  region.  The  arrow  indicates  change  upon  heating 
(22-90  °C). 


(six  for  the  polymer  system  and  four  for  the  gel  and  liquid). 
This  observation  warrants  further  investigation,  e.g.  if  true 
also  for  other  polymers,  gels,  and  liquid  systems. 

However,  mixed  liquid  and  polymeric  environments  of 
a  lithium  cation  are  not  restricted  to  systems  forming  gels. 
Recently,  Wang  and  Balbuena  performed  a  study  of  LiTf  dis¬ 
solved  in  a  mixed  polymer  (PP)  and  liquid  (water)  system 
[12].  The  polymer  acts  as  a  hydrophobic  membrane  allowing 
Li+  transport  to  the  anode,  with  application  in  Li/sea  water 
batteries. 

Mixed  coordination  can  also  result  from  interactions 
with  solids  dissolved  in  polymer  or  liquid  phases,  e.g. 
nano-composite  polymer  electrolytes.  Nano-sized  TiC>2  has 
been  reported  to  enhance  ion  conductivity,  at  least  for 
poorly  conducting  (semi-crystalline)  SPEs,  and  especially 
the  cation  transference  numbers  [113],  while  simultaneously 
increase  the  cyclability.  However,  the  interactions  between 
the  electrolyte  species  and  the  nano-particles  are  not  known. 
Recently,  we  found  a  lithium  cation  to  coordinate  strongly 
to  a  model  rutile  (110)  Ti02  surface  [114].  However,  the 
BF4“  anion  was  the  species  preferably  adsorbed  (Fig.  6),  as 
compared  to  within  a  PEO  matrix.  This  is  in  agreement  with 
NMR  data  showing  slower  anion  dynamics  [115],  and  should 
result  in  increased  cation  transference  numbers,  as  observed. 

As  new  materials  and  species  are  used  as  electrolyte  matri¬ 
ces/solvents  the  lithium  ion  coordination  can  easily  be  mod¬ 
eled  using  ab  initio  calculations  to  assist  in  the  interpretation 
of  the  effects  on  the  ion  conductivity. 

2.4.  Additives 

Organic  solvents  used  in  liquid  electrolytes  are  to  some 
extent  decomposed  during  the  first  charge/discharge  cycle  of 
a  battery  and  a  solid  electrolyte  interface  (SEI)  is  formed. 
Subsequently,  the  SEI  layer  largely  determines  the  perfor¬ 
mance  of  the  electrode/electrolyte.  If  a  small  amount  (a  few 
percent)  of  a  suitable  additive  is  used,  one  which  is  bro¬ 
ken  down  before  the  other  electrolyte  components,  a  better 
SEI  layer  can  be  formed  and  thereby  reducing  further  side- 
reactions  or  exfoliation  of  graphite  electrodes.  Numerous 
experimental  studies  on  both  SEIs  and  additives  exist,  but 
only  a  few  theoretical  studies  are  to  be  found,  most  likely 
due  to  the  more  complex  and  demanding  computational  lev¬ 
els  needed  to  accurately  model  reaction  paths  and  systems 
with  unpaired  electrons  present  (as  for  radical  reactions).  To 
the  best  of  our  knowledge  the  first  study  was  made  by  Scan¬ 
lon  et  al.  [116]  using  ab  initio  calculations  in  an  attempt  to 
correlate  the  stability  of  unsaturated  cyclic  ethers  with  their 
reported  success  to  enhance  cycle  life  by  an  order  of  magni¬ 
tude  in  non-aqueous  electrolytes  [117].  However,  the  authors 
failed  to  find  such  a  relationship,  possibly  due  to  other  factors 
determining  the  experimental  outcome. 

More  recently,  experimental  activity  has  focused  on  differ¬ 
ent  unsaturated  carbonates,  e.g.  vinylene  carbonate  (VC)  and 
vinyl  ethylene  carbonate  (VEC).  The  computational  studies, 
by  DFT  and  using  continuum  methods  to  implicitly  include 
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Fig.  6.  Coordination  of  BF4  to  rutile  (1  1  0). 


a  solvent,  on  VC  and  VEC  has  focused  on  if  the  reduction 
mechanism  is  one-  or  two-electron  based  [118,119],  shown 
that  VEC  has  no  barrier  to  overcome  in  its  path  to  forming 
passivating  lithium  carbonate  (L^COs)  [119],  and  suggested 
how  the  polymerization  of  VC  products  may  take  place,  form¬ 
ing  the  SEI  layer  [118]. 

However,  in  order  to  assess  the  relative  stability  of  the 
additives  versus  the  main  solvents  used,  also  the  decom¬ 
position  of  the  solvents  must  be  modeled.  Endo  et  al. 
modeled  many  different  battery  solvents  using  MNDO  and 
HF  level  calculations  [120,121],  while  later  studies  have 
primarily  studied  EC  and  PC  using  different  DFT  meth¬ 
ods  and  also  included  solvent  effects  via  SCRF  methods 
[90,118,119,122].  We  foresee  this  area  of  applying  ab  initio 
calculations  to  problems  of  battery  electrolytes  to  prosper,  as 
there  are  many  outstanding  issues  with  respect  to  electrolyte 
decomposition  in  general  and  the  formation  and  stability  of 
the  SEI  layer  in  particular. 


3.  Summary  and  conclusions 

The  application  of  ab  initio  calculations  to  various  aspects 
of  battery  electrolytes  and  how  these  can  assist  in  the  design 
of  new  electrolytes  has  been  reviewed.  The  strategy  of  by 
using  ion  pairs  and  anions  to  compute  the  salt  dissociation 
and  the  stability  versus  oxidation  form  a  basis  to  with  con¬ 
fidence  suggest  new  and  better  lithium  salts,  perhaps  based 
on  QSPR  analysis.  The  calculations  can  reveal  the  barriers  of 
conformational  changes  for  polymers  used  to  make  polymer 
electrolytes  and  to  construct  new  potentials  for  MD  simula¬ 
tions.  Furthermore,  the  calculations  can  assist  in  interpreting 
spectroscopic  data,  for  both  liquid,  gel,  and  polymeric  elec¬ 


trolytes,  and  this  way  reveal  the  solvation  of  the  lithium  ion. 
Also  lithium  ion  coordination  within  crystalline  and  nano¬ 
composite  polymer  electrolytes  has  been  modeled.  The  role 
and  use  of  additives  to  liquid  electrolytes  is  an  important 
present  and  future  area.  The  use  of  modeling  techniques  can 
now  follow  almost  any  new  experimentally  launched  concept 
and  contribute  with  accurate  and  cost-effective  data  to  assist 
interpretation. 

Acknowledgements 

The  financial  support  from  the  Swedish  Science  Coun¬ 
cil,  Angpanneforeningens  Forskningsstiftelse,  the  Nordic 
Energy  Research  Council,  the  Swedish  Strategic  Science 
Foundation,  and  Knut  and  Alice  Wallenbergs  stiftelse,  are  all 
gratefully  acknowledged.  The  large  amount  of  computational 
grants  given  to  us  by  the  Swedish  National  Allocations  Com¬ 
mittee  is  gratefully  acknowledged.  The  successful  collabo¬ 
ration  and  many  stimulating  discussions  with  Prof.  Michel 
Armand  have  in  many  aspects  been  crucial  to  our  lithium  salt 
projects. 

References 

[1]  J.-C.  Soetens,  C.  Millot,  B.  Maigret,  J.  Phys.  Chem.  A  102  (1998) 
1055. 

[2]  T.  Li,  P.B.  Balbuena,  J.  Electrochem.  Soc.  146  (1999)  3613. 

[3]  K.  Tasaki,  J.  Electrochem.  Soc.  149  (2002)  A418. 

[4]  M.  Masia,  M.  Probst,  R.  Rey,  J.  Phys.  Chem.  B  108  (2004)  2016. 

[5]  M.  Masia,  R.  Rey,  J.  Phys.  Chem.  B  108  (2004)  17992. 

[6]  PT.  Boinske,  L.  Curtiss,  J.W.  Halley,  B.  Lin,  A.  Sutjianto,  J.  Corn- 
put.  Aided  Mater.  Des.  3  (1996)  385. 

[7]  O.  Borodin,  G.  Smith,  Macromolecules  31  (1998)  8396. 


P.  Johansson,  P.  Jacobsson  /  Journal  of  Power  Sources  153  (2006)  336-344 


343 


[8]  A.  Palma,  A.  Pasquarello,  G.  Ciccotti,  R.  Car,  J.  Chem.  Phys.  108 
(1998)  9933. 

[9]  O.  Borodin,  G.  Smith,  J.  Phys.  Chem.  B  104  (2000)  8017. 

[10]  B.G.  Dixon,  R.S.  Morris,  J.  Power  Sources  119-121  (2003)  856. 

[11]  F.  Miiller-Plathe,  W.F.  van  Gunsteren,  J.  Chem.  Phys.  103  (1995) 
4745. 

[12]  Y.  Wang,  P.B.  Balbuena,  J.  Phys.  Chem.  B  108  (2004)  15694. 

[13]  A.  Hektor,  M.  Klintenberg,  A.  Aabloo,  J.O.  Thomas,  J.  Mater. 
Chem.  13  (2003)  214. 

[14]  M.  Forsyth,  V.A.  Payne,  M.A.  Ratner,  S.W.  de  Leeuw,  D.F.  Shriver, 
Solid  State  Ionics  53-56  (1992)  1011. 

[15]  S.  Neyertz,  D.  Brown,  J.O.  Thomas,  Electrochim.  Acta  40  (1995) 
2063. 

[16]  S.  Neyertz,  J.O.  Thomas,  D.  Brown,  Comput.  Pol.  Sci.  5  (1995) 
107. 

[17]  S.  Neyertz,  D.  Brown,  J.  Chem.  Phys.  104  (1996)  3797. 

[18]  D.  Brandell,  M.  Klintenberg,  A.  Aabloo,  J.O.  Thomas,  Int.  J.  Quan¬ 
tum  Chem.  80  (2000)  799. 

[19]  D.  Brandell,  M.  Klintenberg,  A.  Aabloo,  J.O.  Thomas,  J.  Mater. 
Chem.  12  (2002)  565. 

[20]  C.R.A.  Catlow,  G.E.  Mills,  Electrochim.  Acta  40  (1995)  2057. 

[21]  V.A.  Payne,  J.-H.  Xu,  M.  Forsyth,  M.A.  Ratner,  D.F.  Shriver,  S.W. 
de  Leeuw,  Electrochim.  Acta  40  (1995)  2087. 

[22]  S.  Neyertz,  D.  Brown,  Electrochim.  Acta  43  (1998)  1343. 

[23]  A.  Aabloo,  J.O.  Thomas,  Electrochim.  Acta  43  (1998)  1361. 

[24]  S.W.  de  Leeuw,  A.  van  Zon,  G.J.  Bel,  Electrochim.  Acta  46  (2001) 
1419. 

[25]  A.  van  Zon,  B.  Mos,  P.  Verkerk,  S.W.  de  Leeuw,  Electrochim.  Acta 
46  (2001)  1717. 

[26]  Z.  Gadjourova,  Y.G.  Andreev,  D.P.  Tunstall,  P.G.  Bruce,  Nature 
412  (2001)  520. 

[27]  D.  Brandell,  A.  Liivat,  H.  Kasemagi,  A.  Aabloo,  J.O.  Thomas,  J. 
Mater.  Chem.  15  (2005)  1422. 

[28]  H.  Kasemagi,  M.  Klintenberg,  A.  Aabloo,  J.O.  Thomas,  J.  Mater. 
Chem.  11  (2001)  3191. 

[29]  H.  Kasemagi,  M.  Klintenberg,  A.  Aabloo,  J.O.  Thomas,  Solid  State 
Ionics  147  (2002)  367. 

[30]  H.  Kasemagi,  M.  Klintenberg,  A.  Aabloo,  J.O.  Thomas,  Elec¬ 
trochim.  Acta  48  (2003)  2273. 

[31]  S.D.  Druger,  A.  Nitzan,  M.A.  Ratner,  J.  Chem.  Phys.  79  (1983) 
3133. 

[32]  T.  Matsui,  K.  Takeyama,  Electrochim.  Acta  43  (1998)  1355. 

[33]  J.F.  Snyder,  M.A.  Ratner,  D.F.  Shriver,  J.  Electrochem.  Soc.  148 
(2001)  A858. 

[34]  J.F.  Snyder,  M.A.  Ratner,  D.F.  Shriver,  Solid  State  Ionics  147 
(2002)  249. 

[35]  J.S.  Francisco,  I.H.  Williams,  J.  Phys.  Chem.  94  (1990)  8522. 

[36]  S.P.  Gejji,  K.  Hermansson,  J.  Lindgren,  J.  Phys.  Chem.  97  (1993) 
3712. 

[37]  S.P.  Gejji,  K.  Hermansson,  J.  Lindgren,  J.  Phys.  Chem.  97  (1993) 
6986. 

[38]  S.P.  Gejji,  K.  Hermansson,  J.  Tegenfeldt,  J.  Lindgren,  J.  Phys. 
Chem.  97  (1993)  11402. 

[39]  W.  Huang,  R.  Freeh,  R.A.  Wheeler,  J.  Phys.  Chem.  98  (1994)  100. 

[40]  D.H.  Johnston,  D.F.  Shriver,  Inorg.  Chem.  32  (1993)  1045. 

[41]  B.  Klassen,  R.  Aroca,  G.A.  Nazri,  J.  Phys.  Chem.  100  (1996) 
9334. 

[42]  F.  Ramondo,  L.  Bencivenni,  R.  Caminiti,  C.  Sadun,  Chem.  Phys. 
151  (1991)  179. 

[43]  S.P.  Gejji,  C.H.  Suresh,  L.J.  Bartolotti,  S.R.  Gadre,  J.  Phys.  Chem. 
A  101  (1997)  5678. 

[44]  D.  Benrabah,  R.  Arnaud,  J.-Y.  Sanchez,  Electrochim.  Acta  40 
(1995)  2437. 

[45]  P.  Johansson,  S.P.  Gejji,  J.  Tegenfeldt,  J.  Lindgren,  Electrochim. 
Acta  43  (1998)  1375. 

[46]  I.  Rey,  P.  Johansson,  J.  Lindgren,  J.-C.  Lassegues,  J.  Grondin,  L. 
Servant,  J.  Phys.  Chem.  B  102  (1998)  3249. 


[47]  R.  Arnaud,  D.  Benrabah,  J.-Y.  Sanchez,  J.  Phys.  Chem.  100  (1996) 
10882. 

[48]  S.P.  Gejji,  C.H.  Suresh,  K.  Babu,  S.R.  Gadre,  J.  Phys.  Chem.  A 
103  (1999)  7474. 

[49]  P.  Johansson,  J.  Tegenfeldt,  J.  Lindgren,  J.  Phys.  Chem.  A  104 
(2000)  954. 

[50]  P.  Johansson,  J.  Phys.  Chem.  A  105  (2001)  9258. 

[51]  S.P.  Gejji,  PR.  Agrawal,  N.R.  Dhumal,  Theor.  Chem.  Acc.  107 
(2002)  351. 

[52]  F.  Kita,  H.  Sakata,  A.  Kawakami,  H.  Kamizori,  T.  Sonoda,  H. 
Nagashima,  N.V.  Palenko,  Y.L.  Yagupolskii,  J.  Power  Sources 
97-98  (2001)  581. 

[53]  O.  Borodin,  G.D.  Smith,  R.L.  Jaffe,  J.  Comput.  Chem.  22  (2001) 
641. 

[54]  S.E.  Popov,  A.E.  Nikiforov,  O.E.  Bushkova,  V.M.  Zhukovsky,  J. 
Phys.  Chem.  A  108  (2004)  10280. 

[55]  X.  Zhang,  J.K.  Pugh,  N.  Sukpirom,  M.M.  Lerner,  Int.  J.  Inorg. 
Mater.  2  (2000)  115. 

[56]  C.K.  Kim,  J.  Won,  H.S.  Kim,  Y.S.  Kang,  H.G.  Li,  C.K.  Kim,  J. 
Comput.  Chem.  22  (2001)  827. 

[57]  X.  Xuan,  H.  Zhang,  J.  Wang,  H.  Wang,  J.  Phys.  Chem.  A  108 
(2004)  7513. 

[58]  J.  Grondin,  L.  Ducasse,  J.-L.  Bruneel,  L.  Servant,  J.-C.  Lassegues, 
Solid  State  Ionics  166  (2004)  441. 

[59]  P.  Johansson,  P.  Jacobsson,  J.  Phys.  Chem.  A  105  (2001)  8504. 

[60]  H.  Markusson,  P.  Johansson,  P.  Jacobsson,  Electrochem.  Solid  State 
Lett.  8  (2005)  A215. 

[61]  P.  Johansson,  S.  Beranger,  M.  Armand,  H.  Nilsson,  P.  Jacobsson, 
Solid  State  Ionics  156  (2003)  129. 

[62]  H.  Markusson,  S.  Beranger,  P.  Johansson,  M.  Armand,  P.  Jacobs¬ 
son,  J.  Phys.  Chem.  A  107  (2003)  10177. 

[63]  P.  Johansson,  H.  Markusson,  P.  Jacobsson,  M.  Armand,  Phys. 
Chem.  Chem.  Phys.  6  (2004)  895. 

[64]  P.  Johansson,  in  manuscript. 

[65]  M.  Ue,  A.  Murakami,  S.  Nakamura,  J.  Electrochem.  Soc.  149 
(2002)  1572. 

[66]  J.  Barthel,  H.-J.  Gores,  R.  Neueder,  A.  Schmid,  Pure  Appl.  Chem. 
71  (1999)  1705. 

[67]  D.E.  Fenton,  J.M.  Parker,  P.V.  Wright,  Polymer  14  (1973)  589. 

[68]  F.  Miiller-Plathe,  W.F.  van  Gunsteren,  Macromolecules  27  (1994) 
6040. 

[69]  H.  Yoshida,  I.  Kaneko,  H.  Matsuura,  Y.  Ogawa,  M.  Tasumi,  Chem. 
Phys.  Lett.  196  (1992)  601. 

[70]  M.  Berzaghi,  A.  Gamba,  G.  Morosi,  J.  Mol.  Struct.  (THEOCHEM) 
170  (1988)  69. 

[71]  S.  Tsuzuki,  T.  Uchimaru,  K.  Tanabe,  T.  Hirano,  J.  Phys.  Chem.  97 
(1993)  1346. 

[72]  R.L.  Jaffe,  G.D.  Smith,  D.Y.  Yoon,  J.  Phys.  Chem.  97  (1993) 
12745. 

[73]  G.D.  Smith,  R.L.  Jaffe,  D.Y.  Yoon,  J.  Phys.  Chem.  97  (1993) 
12752. 

[74]  A.  Abe,  T.  Hirano,  T.  Tsuruta,  Macromolecules  12  (1979)  1092. 

[75]  K.  Tasaki,  Polym.  Prepr.  35  (1994)  763. 

[76]  K.  Tasaki,  A.  Abe,  Polym.  J.  17  (1985)  641. 

[77]  A.  Abe,  K.  Tasaki,  J.E.  Mark,  Polym.  J.  17  (1985)  883. 

[78]  G.D.  Smith,  D.Y.  Yoon,  R.L.  Jaffe,  Macromolecules  26  (1993) 
5213. 

[79]  M.A.  Murcko,  R.A.  DiPaola,  J.  Am.  Chem.  Soc.  114  (1992)  10010. 

[80]  F.  Miiller-Plathe,  Acta  Polym.  45  (1994)  259. 

[81]  G.D.  Smith,  D.Y.  Yoon,  R.L.  Jaffe,  Polymer  35  (1995)  680. 

[82]  Y.  Sasanuma,  Macromolecules  28  (1995)  8629. 

[83]  E.E.  Astrup,  Acta  Chem.  Scand.  A33  (1979)  655. 

[84]  S.P.  Gejji,  J.  Tegenfeldt,  J.  Lindgren,  Chem.  Phys.  Lett.  226  (1994) 
427. 

[85]  A.  Sutjianto,  L.A.  Curtiss,  Chem.  Phys.  Lett.  264  (1997)  127. 

[86]  J.-K.  Hyun,  H.  Dong,  C.P.  Rhodes,  R.  Freeh,  R.A.  Wheeler,  J. 
Phys.  Chem.  B  105  (2001)  3329. 


344 


P  Johansson,  P.  Jacobsson  /  Journal  of  Power  Sources  153  (2006)  336-344 


[87]  R.J.  Blint,  J.  Electrochem.  Soc.  142  (1995)  696. 

[88]  R.J.  Blint,  J.  Electrochem.  Soc.  144  (1997)  787. 

[89]  B.  Klassen,  R.  Aroca,  M.  Nazri,  G.A.  Nazri,  J.  Phys.  Chem.  B  102 
(1998)  4795. 

[90]  Y.  Wang,  S.  Nakamura,  M.  Ue,  P.B.  Balbuena,  J.  Am.  Chem.  Soc. 
123  (2001)  11708. 

[91]  P.  Johansson,  M.  Edvardsson,  J.  Adebahr,  P.  Jacobsson,  J.  Phys. 
Chem.  B  107  (2003)  12622. 

[92]  S.-A.  Hyodo,  K.  Okobayashi,  Electrochim.  Acta  34  (1989)  1551. 

[93]  S.-A.  Hyodo,  K.  Okobayashi,  Electrochim.  Acta  34  (1989)  1557. 

[94]  R.  Freeh,  W.  Huang,  Macromolecules  28  (1995)  1246. 

[95]  G.D.  Smith,  R.L.  Jaffe,  H.  Partridge,  J.  Phys.  Chem.  101  (1997) 
1705. 

[96]  D.  Ray,  D.  Feller,  M.B.  More,  E.D.  Glendening,  P.B.  Armentrout, 
J.  Phys.  Chem.  100  (1996)  16116. 

[97]  M.B.  More,  E.D.  Glendening,  D.  Ray,  D.  Feller,  P.B.  Armentrout, 
J.  Phys.  Chem.  100  (1996)  1605. 

[98]  E.D.  Glendening,  D.  Feller,  M.A.  Thomson,  J.  Am.  Chem.  Soc. 
116  (1994)  10657. 

[99]  R.  Freeh,  W.  Huang,  Solid  State  Ionics  72  (1994)  103. 

[100]  R.  Kreitner,  J.  Park,  M.  Xu,  E.M.  Eyring,  S.  Petrucci,  Macro¬ 
molecules  29  (1996)  4722. 

[101]  S.P.  Gejji,  P.  Johansson,  J.  Tegenfeldt,  J.  Lindgren,  Comput.  Pol. 
Sci.  5  (1995)  99. 

[102]  A.  Sutjianto,  L.A.  Curtiss,  J.  Phys.  Chem.  A  102  (1998)  968. 

[103]  P.  Johansson,  J.  Tegenfeldt,  J.  Lindgren,  J.  Phys.  Chem.  A  102 
(1998)  4660. 

[104]  A.G.  Baboul,  PC.  Redfern,  A.  Sutjianto,  L.A.  Curtiss,  J.  Am. 
Chem.  Soc.  121  (1999)  7220. 


[105]  S.P.  Gejji,  S.R.  Gadre,  V.J.  Barge,  Chem.  Phys.  Lett.  344  (2001) 
527. 

[106]  P.  Johansson,  S.P.  Gejji,  J.  Tegenfeldt,  J.  Lindgren,  Solid  State 
Ionics  86-88  (1996)  297. 

[107]  P.  Johansson,  J.  Tegenfeldt,  J.  Lindgren,  Polymer  40  (1999)  4399. 

[108]  P.  Johansson,  Polymer  42  (2001)  4367. 

[109]  PC.  Redfern,  L.A.  Curtiss,  J.  Power  Sources  110  (2002)  401. 

[110]  J.B.  Kerr,  G.  Liu,  L.A.  Curtiss,  PC.  Redfern,  Electrochim.  Acta  48 
(2003)  2305. 

[111]  P.  Johansson,  P.  Jacobsson,  Electrochim.Acta  48  (2003)  2279. 

[112]  J.-M.  Tarascon,  M.  Armand,  Nature  414  (2001)  359. 

[113]  F.  Croce,  R.  Curini,  A.  Martinelli,  L.  Persi,  F.  Ronci,  B.  Scrosati, 
R.  Caminiti,  J.  Phys.  Chem.  B  103  (1999)  10632. 

[114]  P.  Johansson,  P.  Jacobsson,  Solid  State  Ionics  170  (2004)  73. 

[115]  J.  Adebahr,  A.S.  Best,  N.  Byrne,  P.  Jacobsson,  D.R.  MacFarlane, 
M.  Forsyth,  Phys.  Chem.  Chem.  Phys.  5  (2003)  720. 

[116]  L.G.  Scanlon,  W.  Krawiec,  A.T.  Yeates,  D.S.  Dudis,  Electrochim. 
Acta  40  (1995)  2431. 

[117]  K.M.  Abraham,  J.S.  Foos,  J.L.  Goldman,  J.  Electrochem.  Soc.  131 
(1984)  2197. 

[118]  Y.  Wang,  S.  Nakamura,  K.  Tasaki,  P.B.  Balbuena,  J.  Am.  Chem. 
Soc.  124  (2002)  4408. 

[119]  J.M.  Vollmer,  L.A.  Curtiss,  D.R.  Vissers,  K.  Amine,  J.  Electrochem. 
Soc.  151  (2004)  A178. 

[120]  E.  Endo,  M.  Ata,  K.  Tanaka,  K.  Sekai,  J.  Electrochem.  Soc.  145 
(1998)  3757. 

[121]  E.  Endo,  K.  Tanaka,  K.  Sekai,  J.  Electrochem.  Soc.  147  (2000) 
4029. 

[122]  T.  Li,  P.B.  Balbuena,  Chem.  Phys.  Lett.  317  (2000)  421. 


